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The a-stable distributions introduced by Levy play an important role in probabilistic theoretical 
studies and their various applications, e.g., in statistical physics, life sciences, and economics. In 
the present paper we study sequences of long-range dependent random variables whose distributions 
have asymptotic power law decay, and which are called (q, a)-stable distributions. These sequences 
are generalizations of i.i.d. a-stable distributions, and have not been previously studied. Long- 
range dependent (q, a)-stable distributions might arise in the description of anomalous processes 
in nonextensive statistical mechanics, cell biology, finance. The parameter q controls dependence. 
If q = 1 then they are classical i.i.d. with a-stable Levy distributions. In the present paper we 
establish basic properties of (q, a)-stable distributions, and generalize the result of Umarov, Tsallis 
and Steinberg (2008), where the particular case a = 2, q € [1, 3), was considered, to the whole range 
of stability and nonextensivity parameters a £ (0,2] and q £ [1,3), respectively. We also discuss 
possible further extensions of the results that we obtain, and formulate some conjectures. 

PACS numbers: 



I. INTRODUCTION 

The central limit theorem (CLT) and a-stable distri- 
butions have rich applications in various fields including 
the Boltzmann-Gibbs (BG) statistical mechanics. The 
nonextensive statistical mechanics [H-[l| characterized by 
the nonextensivity index q (which recovers the BG theory 
in the case q — 1) studies, in particular, strongly corre- 
lated random states, mathematical models of which can 
be represented by specific long-range dependent random 
variables. The g-central limit theorem consistent with 
nonextensive statistical mechanics was established in pa- 
per Q ■ The main objective of Q was to study the scaling 
limits (attractors) of sums of q- independent random vari- 
ables with a finite (2q — l)-variance. The mapping 



F q -.g q [2\^g z(q) [2] 



(i) 



where F q is the g-Fourier transform (see Section 2) z(s) = 
(1 + s)/(3 — s), and Q q [2] is the set of g-Gaussians up to 
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a constant factor (see, e.g. @, Q), was essentially used 
in the description of attractors. The number 2 in the 
notation will soon become transparent in the context of 
the current paper. 

In the present work we study a g-analog of the te- 
stable Levy distributions. In this sense, the present 
paper is a conceptual continuation of Q. The classic 
theory of a-stable distributions was originated by Paul 
Levy and developed by Levy, Gnedenko, Feller and oth- 
ers; for details and history see, for instance, [8-14] and 
references therein. Distributions with asymptotic power- 
law decay (a-stables, and particularly g-Gaussians, in the 
first place) found a huge number of applications in var- 
ious practical studies (see, e.g., [14-23], just to mention 
a few), confirming the frequent nature of these distri- 
butions. As it will become clear later on, (q, a)-stable 
distributions unify both of them. Indeed, (l,a)-stable 
distributions correspond to the a-stable ones, and the 
(q, 2)-stable distributions correspond to the g-Gaussian 
ones. All (q, a)-stable distributions, except Gaussians 
((1, 2)-distributions), exhibit asymptotic power-laws. In 
practice the researcher is often interested in identification 
of a correct attractor of correlated states, which plays 
a major role in the adequate modeling of physical phe- 
nomenon itself. This motivates the study of sequences of 
(g, a)-stable distributions and their attractors, as focused 



2 



in the present paper. 

For simplicity we will consider only symmetric (q, a 
stable distributions in the one-dimensional case (see 
for the multivariate g-CLT). We denote the class of ran- 
dom variables with (g, a)-stable distributions by C q [a\. 
A random variable X G £ 9 [a] has a symmetric density 

_ l + Q 

f{x) with asymptotics / ~ C\x\ , |a;| — > oo, 

where 1 < g < 2, < a < 2, and C is a positive 
constant. Hereafter g(x) ~ h(x),x — ¥ a, means that 
lim^^a j^pj = 1. Linear combinations and properly scal- 
ing limits of sequences of ^-independent random vari- 
ables with (g, a)-stable distributions are again random 
variables with (q, a)-stable distributions, justifying that 
£ q [a} form a class of "stable" distributions. To this end, 
we note that £ g [a!] shares the same asymptotic behav- 
ior with the set C sym (7) of symmetric Levy distributions 
centered at 0, where 



7 = 7(9, at) 



a(2 - q) 
1 + aO-l)' 



However, there is an essential difference between (q, te- 
stability of g-independent random variables and the 
classic stability of a-stable distributions. Namely, q- 
independence exhibits a special long-range correlation 
between random variables (see the exact definition in 
Section 3). In practice this notion reflects physical states 
(arising e.g. in nonextensive statistical mechanics), which 
are strongly correlated. The term "global correlation" in- 
stead of " strong correlation" is also used widely in physics 
literature; see, e.g. JaQ. Examples of such systems in- 
clude earthquakes [24], cold atoms in optical dissipativc 
lattices [25| , and dusty plasma (2(| . A decomposition of 
nonextensive processes with strong correlation into inde- 
pendent states can not adequately reflect their evolution. 
Likewise, (q, a)-stable distributions can not be captured 
by the existing theory of a-stablc distributions, which is 
heavily based on the concept of independence (or weak 
dependence) . This distinction ends up with essential im- 
plication: attractors of (q, a)-stable distributions are dif- 
ferent from the attractors of a-stable distributions, unless 
q = 1. If q = 1 then correlation disappears, that is q- 
independence becomes usual probabilistic independence, 



and 7(1, a) = a, implying C%[ 



a 



Following the method established in [7(, we will apply 
-Fij-transform for the study of sequences of q-independent 
(q, a)-stable distributions. Parameter q controls correla- 
tion. We will classify (q, a)-stable distributions depend- 
ing on parameters 1 < q < 2 (or equivalently 1 < Q < 3, 
Q = 2q — 1) and < a < 2. We establish the mapping 



F q : g qL [2] -> Q q [a], 



(2) 



where G q [a] is the set of functions {be q ^ , b > 0, f3 > 
0}, and 

l 3 + 

q = -r-; — > <3 = 2<? - 1 , 

1 + a 



1 



l + a(q-l) 



The particular case q = Q = 1 recovers q L = jt^, al- 
ready known in the literature Q. Denote Q\ = {(Q,a) : 
1 < Q < 3, a = 2}, Q 2 = {(Q,a) :1<Q<3, 0<a< 
2} and Q = Qi U Q 2 . Note that the case (Q, a) € Qi for 
g-independent random variables with a finite Q-variance 
was studied in [7|. For (Q, a) £ Q2 the Q- variance is 
infinite. We will focus our analysis namely on the lat- 
ter case. Note that the case a = 2, in the framework of 
this classification like the classic a-stable distributions, 
becomes peculiar. 

In the scope of second classification we study the at- 
tractors of scaled sums, and expand the results of paper 
0] to the region Q generalizing the mapping (J]) to 

: S q [a] -> G Za ( q )[a], 1 < q < 2, < a < 2, (3) 

where 

_ a-2(l-q) _ ag + l-q 

^ayS) — dna z a ys) — — 

a a + 1 — q 

Note that, if a = 2, then (2(q) = q and 22(9) = 
(1 + <z)/(3 — q), thus recovering the mapping (JXJ) , and 
consequently, the result of 0. 

These two classifications of (<?, a)-stable distributions 
based on mappings ^ and ([3]) respectively, can be uni- 
fied to the scheme 



[2] 



(4) 



F a 



which gives the full picture of interrelations of (q, te- 
stable distributions with parameters q £ [1,2) and a € 
(0,2) (see details in Section VIII). 



II. PRELIMINARIES AND AUXILIARY 
RESULTS 

A. Basic operations of g-algebra 

In this section we briefly recall the basic operations of 
g-algebra. Indeed, the analysis we will conduct is entirely 
based on the g-structure of nonextensive statistical me- 
chanics (for more details see [3~5] and references therein) . 
To this end, we recall the well known fact that the classi- 
cal Boltzmann-Gibbs entropy Sbg = ~ 5ZiPi m Pi satis- 
fies the additivity property. Namely, if A and B are two 
independent subsystems, then Sbg {A + B) = Sbg(A) + 
Sbg(B). However, the ^-generalization of the classic en- 
tropy introduced in 1] and given by S q — 1 q ^\ Pi with 
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q E 7Z and Si = Sbg, does not possess this property 
if g 7^ 1. Instead, it satisfies the pseudo-additivity (or 
q-additivity) 0, 0, 01 

S q (A + B) = S g (A) + S q (B) + (l-q) S q (A) S q (B). 

Inherited from the right hand side of this equality, the 
q-sum of two given real numbers, x and y, is defined 
as x ® q y = x + y + (1 — q)xy. The g-sum is commu- 
tative, associative, recovers the usual summing opera- 
tion if q = 1 (i.e. x ©i y — x + y), and preserves 



as the neutral element (i.e. x ©„ 



By inversion, 



we can define the q- subtraction as x Q q y = j-^rrj^ ■ 
The q-product for x, y is defined by the binary relation 



h y 



+ V 



1-9 



Here the symbol [x]+ 



means that [x] + = x if x > 0, and [x] + = if x < 0. 
This operation also commutative, associative, recovers 
the usual product when q = 1, and preserves 1 as the 
unity. The g-product is defined if x x ~ q + y x ~ q > 1. 
Again by inversion, it can be defined the q-division: 



x<Z) q y = (x 



= fx 1 -? 



1)" 



q-generalization of the exponential and cyclic 
functions 



Now let us recall the main properties of two functions, 
g-exponential and g-logarithm, which will be essentially 
used in this paper. Let e q and ln q x denote respectively 
the functions 



= [! + (! — g)x]_j_ q , and ln g x — 



x 1 " 9 - I 



, (x > 0). 



The entropy S q then can be conveniently rewritten in 
the form S q — X^ft hr<? For the g-exponential the 



relations e q qV 



e x e v q and e q +y 



e x (8) g e v q hold true. 



These relations can be written equivalently as follows: 
ln q (x® q y) = ln 9 x+ln g y, and \n q (xy) = (In, x)® q (ln q y). 
The g-exponential and g-logarithm have the asymptotics 



and 



1 + x + |x 2 + o(x 2 ), x ^ 0, 



\n q (l + x) = x - |x 2 +o(x 2 ), x -)• 0, (6) 

respectively. The q-product and g-exponential can be 
extended to complex numbers z — x + iy (see @, HI] ) ■ 
In addition, for q ^ 1 the function e q can be analytically 
extended to the complex plain except the point z — 
— 1/(1 — q) and defined as the principal value along the 
cut (— oo,zo). If q < 1, then, for real y, \e q v \ > 1 and 

\ef\ ~ A 9 (l+y 2 )^>, y ->■ oo, with K q = (l-g) 1 ^ 1 -?). 
Similarly, if q > 1, then < \e q v \ < 1 and \ef \ -J- if 
\y\ — > oo. For complex z it is not hard to verify the power 
series representation 



(5) 



OC 

-,2 



1 + z + ^vA%/, 



z < 



(7) 



where A n (q) = T~[fc=o a k(<l), a k (q) = q - k(l - q). Let 
/, = (— 1/| 1 - q\, 1/| 1 - q\). Then it follows from Q that 
for arbitrary real number x € I q the equation 



n=0 

oo 



(-l)M 2n (g) 
(2n + 2)! ' 



'} 



n=0 



(2n + 3)! 



holds. Define for x £ I q the functions g-cos and g-sin by 
formulas 



1 (x) = l-x 2 J2 



n=0 



{-l) n A 2n (q) 
(2n + 2)! 



(8) 



and 



sin „„., ..-.^H^M^i. ( 9 ) 



£ 

n=0 



(2n + 3)! 



In fact, cos g (x) and sin g (x) can be defined for all real x 
by using appropriate power series expansions. Properties 
of g-sin, g-cos, and corresponding g-hyperbolic functions, 
were studied in (2t|. Here we note that the g-analogs of 
Euler's formulas read 



and 



e q x = cos 9 (x) + i sin g (x), 



e lx + e lx 
cos 9 (x) = , sin g (x) 



21 



It follows from the definitions of cos g (x) and sin g (x), and 
from the equality (e x ) 2 = e 2 f +q y 2 (see Lemma 2.1 in [7]), 



that 



cos ? (2x) = e 2 £l? )x -2sin 2 9 _ 1 (x). (10) 
Denote & q (x) = cos q 2x — 1. Then equation pU|) implies 

= (4-7^ 2 - 1) - 2 sinf^x). (11) 
The following two properties of fy q will be used later on. 

Proposition II. 1 Let q > 1. Then 
1- -2 < $ q (x) < 0; 
2. *g(x) = -2qx 2 + o(x 3 ), x^0. 

Proof. Assume q < 1. Since e 2(J 2 i' < 1, 

then (|11[) immediately implies that ^ q (x) < 0. Fur- 
ther, sin 9 (a;) can be written in the form (see (29l |) 

sin g (x) = p q (x) sm[<p q (x)] : where p q (x) = (e^ 1-9 ^ ) 1/2 
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and (fi q (3 



arctan(l — q)x 



A simple calculation yields 



^q( x ) > — 2. Using the asymptotic relation ([5]), we get 

efq-i )x2 - 1 = 2(1 - q)x 2 + o(x 3 ), x^O. (12) 

In turn, it follows from ([9]) that 

- 2 sinf g _ x (a;) = -2 x 2 + o(x 3 ), x -> 0. (13) 

Now dTTJ, (JTSJl and imply the second part of the 
statement. ■ 

Representation shows the behaviour of g- 
cxponential near the origin. It is not hard to verify that 
in the case q > 1 for x > (q — 1) _1 the representation 



r* = [( ? - i)x]-*=i{i - 



i 



(1 - qfx 



1 °° 
fl - qYx 2 53 7^ 



(-l) n ^n(?) 



(1 - g) 4 jj 2 ^ (n + 2)!(g - l) 2n y x 



holds. 



C. g-Fourier transform for symmetric densities 

The q-Fourier transform for q > 1 was introduced in 
and used as a basic tool in establishing the g-analog 
of the standard central limit theorem. Formally the q- 
Fourier transform for a given function f(x) is defined by 



F q [fK0 



ei x( ® q f{x)dx . 



(14) 



For discrete functions fk,k = 0,±1,..., this definition 
takes the form 



(15) 



k— — c 



In the future we use the same notation in both cases. 
We also call (fl"4"|) or (|15p the q- characteristic function of 
a given random variable X with an associated density 
f(x), using the notations or F q [f] equivalently. 

It should be noted that, if in the formal definition (|T4|) 
/ is compactly supported then integration has to be taken 
over this support, although, in contrast with the usual 
analysis, the function e l q ^® q f{x) under the integral does 
not vanish outside the support of /. This is an effect of 
the q- product. 

The g-Fourier transform for nonnegative f{x) can be 
written in the form 



F q [f}(0 = 



dx. 



(16) 



We note that, if the g-Fourier transform of f(x) defined 
by (fl"4"|) exists, then it coincides with (fT5|) . The q- Fourier 
transform determined by the formula (|16[) has an advan- 
tage to compare to the formal definition: it does not use 
the (/-product, which is, as we noticed above, restrictive 
in use (for q > 1). 



Proposition II. 2 Let f(x) be an even function. Then 
its q-Fourier transform can be written in the form 



F g [fK0 



f(x)cos q (x^f(x)]^ i )dx. (17) 



Proof. Notice that, because of the symmetry of /, 



® q f(x)dx . 



e q xi ® g f{x)dx 
unt, ' 

4* ®, /(a) + e" te * ® q f{x)) dx. 



Taking this into account, we have 
1 f°° 

F q [f](Z) = 2j_ a 
Now due to (fTBT) we obtain 



F q [f}(0= /(*)- 



ixZ\f(x)Y 



-ix£[/(x)]« 



-dx . 



which coincides with (|17p. 
Further, denote 



Hn 



{/ g Li : f(x) ~ C\x\-^=^, \x\ -> oo}. 



For a given / G Ff qyCC the constant C = Cj is defined 
uniquely by /. It is readily seen that 0(g, a) 



> 



l+a(q-l) 

1 for all a G (0, 2) and g G [1,2). Moreover, 0(g, a)(2g - 
1) < 3 for all a £ (0,2) and q G [1,2), which implies 
a 2 q -i(f) = 00 • Notice also <fi(q,a) = 1 + a*(g,a), where 



a* = a*(g, a) 



a(2-g) 
l + a(g-l) 



On the other hand, the density <? of any a*-stable 
Levy distributions has the asymptotic behaviour g(x) ~ 
C/\x\ 1+a * , \x\ -> oo. Hence, for a fixed g G [1,2) iIg, Q 
is asymptotically equivalent to the set of densities of a*- 
stable Levy distributions. 

The following proposition plays a key role in our fur- 
ther analysis. 

Proposition II. 3 Let f{x), x G R, be a symmetric prob- 
ability density function. Further, let either 

(i) the (2g — \)-variance cr|g-l(/) < 00 j (associated 
with a — 2, and 1 < g < 2j, or 

(iij /(#) € H q . a , where [2q— l,a) G Q.%. 

Then, for the q-Fourier transform of f(x), the follow- 
ing asymptotic relation holds true: 

F q [fm = i - n 9 , a \t\ a + oaeru -> o, (is) 

where 

if a = 2 ; 



2 2 -°(l+ a (g-l))C / 



Jo 



2-g JO 2/° 

^^-l(/) = /" oc [/W] 29 - 1 ^ 



dy, tf(2«/-l,a)e Q 2 
(19) 
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Proof. First, assume that a = 2. By Proposition 111.21 

/>oo 

F q [fKO= I (e^)® q f(x)dx 

f(x)cos q (xaf(x)] q - 1 )dx. (20) 



Making use of the asymptotic expansion ([5]) we can 
rewrite the right hand side of (|2U)) in the form 

/oo 
mii+ixziKx)]*- 1 
-oo 

- q /2x 2 e[m] 2{q - l) }dx+o{e) 

= 1 - ((7/2)eV 2 Vi^-i + o(C 3 ), £ 0, (21) 

from which the first part of Proposition follows. 

Now, we assume (2q— 1, a) G Qi- We apply Proposition 
IL2] to obtain 

/OO 
/(zMcos^/^r 1 )-!]^ 
-OO 

+ 2f/(,)f,(^)&, 



where AT is a sufhciently large finite number. In the first 
integral we use the asymptotic relation 



o(x 3 ), a; — > 0, which follows from Proposition III. 11 and 
E \j N f(x)%C-^l)dx = 

-q? x 2 f- 1 (x)dx + o(e),^0, (22) 
Jo 

that is a quantity of order o( |^| 5 ) , ^ — > 0, for any <5 < 2. 
In the second integral taking into account the hypothesis 
of the proposition with respect to f(x), we have 



Jn z 



(23) 



= 2C t 



(a + l)(,-l) 



We use the substitution 



00 1 i+^-i) £ 



2?y 



in the last integral, and obtain 
Jn z 



where 



» q ,a\ti\ a + o(\Zn, €->0, (25) 



2 2 -«(l + a( 9 -l))C / /■» ¥„(!,) 
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o V 



Q + l 3 



Hence, the obtained asymptotic relations (|22|) (we take 
<5 G (a, 2)) and (|2"5)l complete the proof. ■ 

For stable distributions /i gjQ , must be positive. We have 
seen (Proposition lll.lj ) that if q > 1, then 4 , 9 (x) < 
(not being identically zero), which yields /i 9jQ > . Note 
also that the condition for f(x) to be symmetric was not 
required in Q, if <T2 9 -i(/) < °°- 



III. WEAK CONVERGENCE OF CORRELATED 
RANDOM VARIABLES 

Let us start this section by introducing the notion of 
q-independence. We will also introduce two types of con- 
vergence, namely, q-convergence and weak q-convergence 
and establish their equivalence. 

By definition, two random variables X and Y are said 
to be (q ,q,q ) -independent if 



F q , [X + Y](0 = F q [X}(0 J^PKO 



(26) 



In terms of densities, equation (|26|) can be rewritten as 
follows. Let /x an d fy be densities of X and Y respec- 
tively, and let fx+Y be the density of X + Y. Then 

V 5 %> fx +Y (x)dx = F q [f x ](£) ® q// F q [f Y U). 

_ J 

(27) 

If all three parameters q , q and g coincide, i.e. q = q = 
q , then we call simply q- independent. For q = 1 the 
condition (12611 turns into the well known relation 



F[fx * fy] = F[f x ] ■ F[fy] 

between the convolution (noted *) of two densities and 
the multiplication of their (classical) characteristic func- 
tions, and holds for independent X and Y. If q ^ 1, 
then (q ,q,q )-independence describes a specific class of 
correlations. 

Remark III.l It is worth to mention at this point that 
g-independence appears to be relevant to the notion of 
scale-invariance [391 ]. To be more specific, it might 
well be that q-independence implies scale-invariance, 
i.e., scale-invariance is necessary for q- independence, al- 
though it is by now clear that it is not sufficient. Indeed, 
scale-invariant probabilistic models exist in the litera- 
ture. Some of them presumably involve q-independence 
since their N — > oo limits are g-Gaussians [39|, |4fJ; oth- 
ers do not involve ^-independence [4ll - |43T | (if they did in- 
volve, their N — » oo limits would have to be g-Gaussians 
and they are not). See also paper [44| which discusses 
limit distributions in general setting within the exchange- 
ability concept. 

Let X/v be a sequence of identically distributed random 
variables. Denote Yn = Xi + ...+Xn. By definition, Xn 
is said to be (q ,q,q ) -independent (or (q ,q,q )-i.i.d.) if 
the relations 



F,[Y N ](0=F q [X 1 ](0 



F g [X N }(0 (28) 
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hold for all TV = 2, 3, .... 

For q = q = q =1 the condition (|28p turns into 
the condition for the sequence X^ to be usual i.i.d. If 
q = q = q then we call the sequence Xn simply a q- 
i.i.d. Consider example of an (q , q, g)-i.i.d. sequence of 
random variables, where q G (1, 3) and q = (3g — l)/(q + 
1). Assume Xjy is the sequence of identically distributed 
random variables with the associated Gaussian density 

G g <(/3,*) = ^e/* 2 , 

where C q i is the normalizing constant (see, e.g. 0]). 
Further, assume the sums X\ + ... + X/v, N — 2,3, 

are distributed according to the density Gj (a, x), where 

i 

a = N 2 -i' j3. Then the sequence X N satisfies l[2"8"|) 
for all N — 2,3,..., with q = q , thus being (q ,q,q)- 
independent identically distributed sequence of random 
variables. 

For the sake of simplicity in this paper we will consider 
only g-i.i.d. random variables. 

By definition, a sequence of random variables Xn is 
said to be q-convergent to a random variable Xqo if 
limN^oo F q [X N }(0 = F q \Xoa\(£) locally uniformly in £. 

Evidently, this definition is equivalent to the weak con- 
vergence (denoted by "=>") of random variables if q = 1. 
For q 7^ 1 denote by W q the set of continuous functions 

satisfying the condition \(f>(x)\ < C(l + i- 1 , x G R. 

A sequence of random variables Xn with the density 
/jv is called weakly q-convergent to a random variable 
Xoo with the density / if J R fN(x)dm q J Rd f(x)dm q 
for arbitrary measure m q defined as dm q {x) — <j) q {x)dx, 
where 4> q G W q . We denote the g-convergence by the 
symbol =l>. 



ir \ R J-R 



- I {l-e\ xU * )di\d l x N {x). (29) 



It is not hard to verify that 



1 



R 



2 sinj_ (Rx(2 - q)t) 



Rj- R q K Rx(2-q)t 



It follows from (J29]) and ([30]) that 

r-R 



(30) 



^ J B (l - F q [f N }(0)dC 



- sin^(x(2 - q)Rf q N 1 Y 



1 - 



Rx(2 - q)^ 1 



d/j. N {x). (31) 



Since 1 < q < 2 by assumption, — > 1 as well. 

It is known [29l - l3~i| that for any q > 1 the prop- 
erties sin q i(x) < 1 and (sin q >(x))/x — > 1, x — > 
hold. Moreover, (sin^ (a;))/x < l,Vx G i?. Suppose, 
Hm^i^oolxl/^ -1 = L N: N > 1. Divide the set {TV > N Q } 
into two subsets A = {Nj > Nq : Ljy > 1} and 
B = {N k > N : L N . < 1}. If N G A, since sin^L. < 1, 

2-q 

there is a number a > such that 



> 2 



X |> V i?|a:|(2 -<?)/« 



— djjL N {x) 



N 



Proposition III. 2 Let q > 1. TTien JV^ 
4> Xq. 



Xq yields 



> Cfx N (\x\ > a) , C > ViV G A, 

for i? small enough. Now taking into account the q- 
convergence of Xm to X and, if necessary, taking R 
smaller, for any e > 0, we obtain 



1 r R 

The proof of this statement immediately follows from ^ N > a ) < / (l _ F q [f ](£))d£ < e, ViV G A 

e obvious fact that W„ is a subset of the set of bounded CR J_ R 



the obvious fact that W q is a subset of the set of bounded 
continuous functions. Recall that a sequence of prob- 
ability measures /Ujv is called iig/ii if, for an arbitrary 
e > 0, there is a compact X e and an integer N* such 
that fi^(R d \ K e ) < e for all N > iV*. 

Proposition III. 3 Lef 1 < g < 2. Assume a sequence 
of random variables Xn , defined on a probability space 
with a probability measure P, and associated densities 
/at, is q-convergent to a random variable X with an as- 
sociated density f. Then the sequence of associated prob- 
ability measures /ijv = PiX^ 1 ) is tight. 

Proof. Assume that 1 < q < 2 and Xn is a q- 
convergent sequence of random variables with associated 
densities /jv and associated probability measures ^n- We 
have 



If N G B, then there exist constants b > 0, S > 0, such 
that 

M*)<TT<^r. \x\>b,VN€B. 
\x\i- 1 \x\i~ 1 

Hence, we have 



IJ,n(\x\ > b) = I f N (x)dx 

\x\>b 



<(1 + S) 



dx 



— , N e B. 



\x\>b \x\ i- 1 



LJ\i-F q [f N m)dt=±J\i- 



fNe q 



dx)d£ 



Since, l/(q—l) > 1, for any e > we can select a number 
b e > b such that /i/v(M > b £ ) < e, N e B. As far as A U 
B = {N > Nq} the proof of the statement is complete. 
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Further, we introduce the function 

D q {t) = D q (t; a) = ief ^ 

= t(l + i(l-q)at q - 1 y^, (32) 

defined on [0, 1], where 1 < q < 2 and a is a fixed real 
number. Obviously, D q (t) is continuous on [0, 1] and dif- 
ferentiable in the interval (0, 1). In accordance with the 
classical Lagrange average theorem for any t 1 ,t 2 , < 
ti < t 2 < 1 there exists a number t± < i* < t 2 such 
that 

D q (h) - D q (t 2 ) = D' q (U)(ti - t 2 ), (33) 

where D q means the derivative of D q {t) with respect to 
t. 

Consider the following Cauchy problem for the 
Bernoulli equation 

1 ia(q — 1) 



V -p 



f 



-y\ i/(0) = o, 



(34) 



It is not hard to verify that y(t) = D q (t) is a solution to 
problem (|34|) . 



Proposition III. 4 For D q (t) the estimate 

\D q (t;a)\<C(l + \a\)-^, t e (0, 1], a e R 1 , (35) 
holds, where constant C does not depend on t. 
Proof. It follows from ([32]) and ((34]) that 
\y'(t)\ <t- l \y + ia{q-l)y«\ 

= |ef _1 +ia( 9 -l)t»- 1 (ej ,t, " 1 )«| 
= |1 + io(l - g)i 9_1 |~^ 

<C(i + |a|)-A,te(o,i]. 

■ 

Now we are in a position to formulate the following 
two theorems on the relationship between ^-convergence 
and weak g-convergence. 

Theorem III. 5 Let 1 < q < 2 and a sequence of random 
vectors Xn be weakly q-convergent to a random vector X . 
Then Xn is q-convergent to X . 

Proof. Assume Xn, with associated densities /at, is 
weakly q-convergent to a X, with an associated density 
/. The difference F q [fN]{Q — J~q[f](0 can be written in 
the form 



(D q (f N (x))-D q (f(x)))dx, (36) 



where D q (t) = D q (t; a) is defined in (|32p with a = x£. It 
follows from (|3"3"|) and ([33)) that 

l-^[/iv](0-^[/](OI 



<C/ 1(1 + (/tv(^) - /(x)) 

which yields F q [f N ](£) ->■ ^[/jv](0 for all £ G I 



Theorem III. 6 Lei 1 < q < 2 and a sequence of ran- 
dom vectors Xn with the associated densities /at is q- 
convergent to a random vector X with the associated den- 
sity f and .Fg [/](£) is continuous at £ = 0. Then Xn 
weakly q-converges to X . 

Proof. Suppose that /jy converges to / in the sense 
of g-convergence. It follows from Proposition IIII.3I that 
the corresponding sequence of induced probability mea- 
sures hn = P(X N 1 ) is tight. This yields relatively 
weak compactness of hn- Theorem lHI . 5 1 implies that each 
weakly convergent subsequence {^Nj} of /Mr converges to 
(i = P{X^ 1 ). Hence, /ijy => £«, or the same, Xn => X. 
Now applying Proposition lIII.2l we complete the proof. ■ 



IV. SYMMETRIC (q, a)-STABLE 
DISTRIBUTIONS AND THEIR PROPERTIES 



In this section we introduce the symmetric (g, a)-stable 
distributions and classify them on the base of mapping 
@. In this classification q takes any value in [1,2), how- 
ever we distinguish the cases a — 2 and < a < 2. 



Definition IV. 1 A random variable X is said to have 
a (g, a) -stable distribution if its q-Fourier transform is 

represented in the form e q /3 '^' , with f3 > 0. We denote 
the set of random variables with (g, a) -stable distributions 
by C q [a]. 

Denote G q [a] = {be q /3 ' 5 ' , b > 0, (3 > 0}. In other words 
X e C q [a], if F g [X] S Gg[a] with 6 = 1. Note that 
if a = 2, then Gq[2] represents the set of g-Gaussians 
and C q [2) - the set of random variables whose densities 
are g*-Gaussians, where g* = (3g — 1)/(1 + g). Fur- 
ther, from the asymptotic relation ([5]) we have e q = 
1 - P\C\ a + o(|C| a ). This and Proposition QX3] imply that 
the associated density of any (g, a)-stable distribution 
belongs to H q ^ a . 

Proposition IV. 2 Let q-independent random variables 
Xj e £q[a],j = 1, .., m. Then for constants 

m 

J2 a i x 3 <= C q [a\. 

2=1 

Proof. Let 

F q [X j }($=e q t\Z\ a ,j = l,...,m. 

Using the properties e q ® q e v q = e q +v and F q [aX](£) = 
F q [X](a 2 ~ q t;), it follows from the definition of the q- 
independence that 

m m 
2=1 2=1 
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Proposition IIV.2I justifies the stability of distributions 
in C q [a] . Recall that if q — 1 then g- independent random 
variables are independent in the usual sense. Thus, if 
q = 1, < a < 2, then Ci[a) = C sym [a], where £ sym [a] 
is the set of a-stable Levy distributions. 

Moreover, the appropriately scaling limit of sequences 
of ^-independent random variables with (g, a)-stable dis- 
tributions has again a (q, a)-stable distribution. To this 
end consider the sum 



1 



J N 



s N (q,a 



X N ),N = 1,2, 



where sj\r(g,a) is a scaling parameter specified below. 
First we prove a general result. 

Theorem IV.3 Assume (2g — 1, a) £ Q 2 . Let X N be 

symmetric q -independent random variables all having the 
same probability density function f(x) £ H qa . Then 

Zn, with SN(q,ce) = (fi q . a N) a ( 2 -9) , is q-convergent to 
a (q, a)-stable distribution, as N — > oo. 

Proof. Assume {Q,a) £ <2 2 - Let / be the density as- 
sociated with X%. First we evaluate F q [Xi] = F q [f(x)]. 
Using Proposition III. 31 we have 



F g [f}(0 = 1 - ^,«|£| Q + °(l£| a U -> 0. (37) 



Denote K 



N~ 



Xj,j = 1,2,.... Then Z N = Y x + 
... + Yn ■ Further, it is readily seen that for a given ran- 
dom variable X and real a > 0, the equality F q [aX](£) = 
F q [X](a 2 ~ q t;) holds. It follows from this relation that 
F q [Y 3 ] = F q [f]{ { J N)1/a ), j = 1,2,... Moreover, it fol- 
lows from the g-independence of X\,X 2 , ■ and the as- 
sociativity of the g-product that 



F q [Z N }(0 
= F q [f]( 



{fJ, q>a N)c . ' [Vq,otN) 

n factors 



-). (38) 



Further, making use of the expansion 
logarithm, equation (|38j) implies 



for the q- 



\n q F q [Z N }(0 = Nln q F q [f}((^ a N)--0 

|£|a ICIa 

= -|£| a + o(l),i\T->oo, (39) 
locally uniformly by £. Hence, locally uniformly by £, 



JV 



lim F q [Z N ] = e"l«l a eG q [a]. 



(40) 



Thus, Zjv is g-convergent to a random variable with 
(g, a)-stable distribution, as N — > oo. ■ 

Since the density of X £ C q [a] is in H q [a] it follows 
immediately the following Corollary from Theorem IIV. 31 



Corollary IV. 4 Assume (2q — 1, a) £ Q 2 - Let Xn be 

a sequence of symmetric q-independent (q, a) -stable ran- 
dom variables. Then Zn, with the same SAr(g,a) in The- 
orem UVlh q-weakly converges to a (q,a)-stable distribu- 
tion. 

Note that a — 2 is not included to Q 2 in Theorem 
IIV. 31 The case a = 2, in accordance with the first part of 
Proposition 111.31 coincides with Theorem 2 of pj. Recall 
that in this case C q [2] consists of random variables whose 
densities are in Q q * [2] , where q* = . 

Theorem II V. 3 1 also allows to establish a connection be- 
tween the classic Levy distributions and g^-Gaussians. 
Indeed, for a X £ £ q [a], its density function / has 
asymptotics 



/ ~ C f /x 



(a+l)/(l+a( 9 -l)) 



\x\ — ► OO. 



It is not hard to verify that there exists a g^-Gaussian, 
which is asymptotically equivalent to /. Let us now find 
Any g^-Gaussian behaves asymptotically Ci/|a;|' ) = 



CVM 2 /^- 1 ), Cj = const, j-. 
Hence, we obtain the relation 

a + 1 



1,2, i.e. 77 = 2/(^-1) 



(41) 



l + a(q-l) qt-l 
Solving this equation with respect to g„ , we have 
L _ 3 + Qa 



Q a = 



a + 1 



-, Q = 2q-1, 



(42) 



linking three parameters: a, the parameter of the te- 
stable Levy distributions, g, the parameter of correla- 
tion, and q%, the parameter of attractors in terms of 
Gaussians. Equation (|42|) identifies all (Q, a)-stable dis- 
tributions with the same index of attractor G q L , proving 
the following proposition. 

Proposition IV.5 Let I < Q < 3 (Q = 2q- 1), < a < 
2. and 



3 + Qa 
a + 1 



Qa . 



(43) 



Then the density of X £ C q [a] is asymptotically equiva- 
lent to q%- Gaussian. 

In the particular case Q = 1, we recover the known 
connection between the classical Levy distributions (g = 
Q = 1) and corresponding g^-Gaussians. In fact, putting 
Q = 1 in equation (|4"2")l . we obtain 



1 



< a < 2. 



(44) 



When a increases between and 2 (i.e. < a < 2), g^ 
decreases between 3 and 5/3 (i.e. 5/3 < g^ < 3). 

It is useful to find the relationship between n — i_ x , 
which corresponds to the asymptotic behaviour of the 
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attractor depending on (a,Q). Using formula (|4"Tj) . we 
obtain 



V = 



2(« + l) 
2 + a(Q-l)' 



(45) 



Proposition IV.6 Let X e C Q [a], 1 < Q < 3, < a < 
2. Then the associated density function fx has asymp- 



totics fx{x) ■ 
fined in J^5[ ). 



\x\ — ¥ oo, where rj = rj(Q, a) is de- 



li Q = 1 (classic Levy distributions) , then (|45|) implies 
the well-known fact rj = a + 1 . 

Analogous relationships can be obtained for other val- 
ues of Q. We call, for convenience, a (Q, a)-stable dis- 
tribution a Q-Cauchy distribution, if a = 1. We obtain 
the classic Cauchy-Poisson distribution if Q = 1. For 
Q-Cauchy distributions (|4"3"|) and (|4"5)l imply 



9i (Q) = -^r - and ^ = 



Q + r 



(46) 



respectively. 



V. SCALING LIMITS OF SUMS OF 
(q, a)-STABLE DISTRIBUTIONS 

In this section we generalize the g-central limit theorem 
established in 7] for g-Gaussians, that is in the case of 
a = 2, to symmetrical (g, a)-stables with any a S (0, 2]. 

Let 1 < q < 2, and / G Q q [a], < a < 2. It fol- 
lows from the definition of the g-exponential that / ~ 
Cf {xl^ 71 , Cf > 0, as |x| — > oo. Analogously, if g <E Q q [2], 

-2 

then g ~ C g \x\ i- 1 , C g > 0, as \x\ — > oo. Comparing or- 
ders of asymptotics we can easily verify that for a fixed 
a € (0, 2] and for any q g (1, 2) there exists a one-to-one 
mapping 

M q , q . : S 9 [a] -+ Q q , [2] , g* = a + 2(g ~ 1} , 

a 

such that the image of a density / e is again den- 

sity. Analogously, there is a one-to-one mapping 

with the same g*, such that it maps /(x) = e q ^ x ^ , 
an element of Q q [a] with the coefficient 6=1 onto the 

element g(x) = e q , 2 with the same coefficient 6=1. 
We notice that if a = 2, then g* = g and both operators 
coincide with the identity operator. 



Let J- q be an operator defined 



To 



fC z(q'),q, F Q' M 'i<i^ wnere z (l*) = irfr- K is read- 
ily seen that in the particular case a = 2 it coincides 
with the g-Fourier transform, T q 
generalized g-Fourier transform. 



F q . We call T q a 



Proposition V.l Assume < a < 2 and Zei i/ie num- 
6ers q* , g» and g 6e connected through the relationships 



. a-%-1) 
g = ana g» 



Qg+ (g- 1) 
a + (°-9)- 



Then the mapping 



T q ■ G q [a] -t Q q , 



(47) 



(48) 



holds. 



Proof. We use the scheme 



Q q [a] Qq, [a] 



(49) 



g q ,[2] -Ag z ^[2] 



for the proof. Let a density / G <?q [a], i.e. asymptotically 
/(x) ~ C/ l^l-a/C'?-!)^ x oo with some C f > 0. Its 
image -Mq,q*[f](x), a g*-Gaussian G q *(f3;x), in order to 
be asymptotically equivalent to /, necessarily 



G q .{[3 V x) 



Ci 



Cr 



x « 



x — > oo. 



Hence, 



a + 2(g-l) 



= 1 



2(9-1) 



a a 
Further, it follows from Corollary 2.10 of Q, that 

F q . : £,,[2] ->a 9l [2], 

where 

1 + g* a+(g-l) 



'7i 



3 — g* a — (g — 1) 



Now taking into account the asymptotic equality (the 
right vertical line in (14910 



G qi (Pi;x) 

we obtain 
9 



C 2 



C 3 



x "I" 1 



!<«-' 



00, 



aq~(q- 1) 



a - (g- 1) 



«(g - 1) 



Thus, the mapping (|48l) holds with g* and g* in equation 
63). ■ 

Let us now introduce two functions that are important 
for our further analysis: 



as — (s — 1) a(s — 1) 

z a\S) — -, —- = 1 H -, 7T : 

a — (s — 1) a — (s — 1) 



(50) 
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where < a < 2, s < a+1, and 

* + 2(s-l) 2(s-l) 



Us) 



1 



a 



, < a < 2. (51) 



It can be easily verified that £ Q (s) = s if a = 2. 

The inverse, z~ x (t), t € (1 — a,oo), of the the first 
function reads 



at+(t- 1) 
a + {t-l) 



1 



oj(t - 1) 
*+ (i- 1) 



(52) 



The function z(s) possess the properties: z a ( - ) 
i and z Q (i) = 



If we denote q Qi i = z Q (g) and 



q a , 



^a 1 ^), then 



z a (— ) = - and z (-) 
9a, l 9 9 



(53) 



Proposition IV. II implies that for < a < 2 and 1 < 
q < min{2, 1 + a} the following mappings hold: 

(i) T 3 : </ g [a], 

(ii) •T 7 " 1 : Zc ,( 9 )N -> £/ 9 [a], 

where J 7 ^ 1 is the inverse to T q . 

It should be noted that as Hilhorst [IH noticed q- 
Fourier transform in general is not one-to-one in the space 
of densities. In paper [33[ the invertibility of F q in the set 
of g-Gaussians is established. Since mappings M. q>q * and 
IC qiq * are one-to-one, relationship (|49p yields invertibility 
of T q in Gz a (q)[ a ] an d validity of property (ii). 

Further, we introduce the sequence q a , n — z a , n {q) = 
z(z a ,n-i (?)),« = 1,2,..., with a given q = z (q), q < 
1 + a. We can extend the sequence g Qj „ for negative in- 
tegers n = —1, —2, ... as well, setting q a ,-n = z a ,~ n {q) = 
z~ (z a ,i- n (q)), n — lj 2, ... . It is not hard to verify that 



1 + 



a(q — 1) aq — n(q — 1) 



a — n(q — 1) a — n(q — 1) 



(54) 



for all integer n satisfying — oo < n < The re- 

striction n < [a/(q — 1)] implies the necessary condi- 
tion q a ^ n > 1, since q- Fourier transform is defined for 
q > 1. Note that q Qj „ is a function of (q,n/a), that 
q a n = 1 for all n = 0,±1,±2, if q = 1, and that 
lim„_>.± oo z a ^ n (q) — 1 for all q ^ 1. Equation (|54j) can be 
rewritten as follows: 



qa.,n 1 



g-i 



0,±1,±2, 



(55) 



We note that the latter coincides with equation (13) 
of (HI, once we identify a with the quantity z therein 
defined, which was obtained through a quite different 
approach (related to the renormalization of the index q 
emerging from summing a specific expression over one 
degree of freedom). 

We also note an interesting property of q a ^ n . If we 
have a (/-Gaussian in the variable |x| Q / 2 (q > 1), i.e., a q- 
exponential in the variable \x\ a , its successive derivatives 



and integrations with respect to |a;| a precisely correspond 
t° <? Q .ri-exponentials in the same variable \x\ a . 

Further, we introduce the sequence q* an — ((q a ,n), 
which can be written in the form 

* -la ^g" 1 ) - a+(n-2)(l-g) 
q «>"- L+ a-n(q-l) ~ a-n(q- 1) 1 [bb) 

for n = 0, ±1, or, equivalently, 

2 a 

n= r, n = 0,±l,..„ (57) 



q* a ,n - 1 q - 1 



It follows from Proposition IV. II and definitions of se- 
quences q a . n and g* n that 

F qa ,n ■ G qa Ja] -> fc|B+1 ,-oo < n < [— []• (58) 



Proposition V.2 For aH n = 0, ±1, ±2, ... i/ie following 
relations 



1 



9a,n+l 



92,n = 52,n 7 

hold. 

Proof. We notice that 



(59) 



(60) 



1 



= 1 



2(9-1) 



9 Q ,n+l a -(rc- l)(g-l) 

On the other hand, by (151)1) 
2(9-1) 



a — (n — l)(g — 1) 



-i * 

^ 9a, ra — 1 1 



which implies (|5"9"|) immediately. The relation (IStJl) can be 
checked easily. ■ 

The property n = 92, n shows that the sequences (|54|) 
and ([55)) coincide if a = 2. Hence, the mapping (|58")) 
takes the form Kj, : £g 2 ,n(2) — > ^g 3 ,„+i(2) ; recovering 
Lemma 2.16 of [7|- Moreover, in this case the duality 
(|59|) holds for the sequence g Q: „ as well. If a < 2 then 
the values of q* an are distinct from the values of q a , n - 
The difference is given by 



q, 



a,n Ha.n 



(2 -«)(!- 9) 
a + n(l — q) 



vanishing for a — 2 , Vg, or for 9 = 1 , Va. In the latter 
case q atn = q* a n = 1. 

Further, we define for n = 0,±1,..., k = 1,2,..., n + 
k < + 1, the operators 
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and 

Fn k (f)=F- a 1 n _ k 0...oT- a 1 n _ 1 [f] 

= ^ 1 „-J-^1-J^1- 1 [/]]-]■ 

In addition, we assume that Fi\f \ = f, if k = for any 
appropriate q. Summarizing the above mentioned rela- 
tionships, we obtain the following assertions. 

Proposition V.3 The following mappings hold: 

1- Fq a>n ■ Qq a , n [a] -> M, -oo < n < [^-]; 

«■ Fn-G^.M^G^M k = 1,2, 

n = 0,±1, ...,-oo < n + k < [^-] + 1; 

3. limk^-oo ^g q [a] = g[a], n = Q,±l,..., 

where G{a) is the set of densities of classic symmetric 
a-stable Levy distributions. 



use of the properties of the q-logarithm, from (|6Tj) we 
obtain 

1*W F qa k [Z N m = N\n qa k F qa k [f}(- 1— j-) 

= N\n Q .(1- l£i- +0 (l£L 
= -|^| Q +o(l), JV->oo, (62) 
locally uniformly by £. Consequently, locally uniformly 

i'.v c 

lirn^ F q ^ k (Z N ) = e-Jf g g qa _ k (a) . (63) 

Thus, Zpf is g Q ,A;-convergent. 

To show the second part of Theorem we use Proposi- 
tion IV31 In accordance with this lemma there exists a 

_|A|Q 

density f(x) G Gq clik _ 1 [a], such that J r a k _ 1 [f] = e g J^ . 
Hence, Zn is g^-convergent to a (g a) jt_i, a:)-stable dis- 
tribution, as N —¥ oo. ■ 



Theorem V.4 Assume < a < 2 and a sequence 
Qa,n, ~°o < n < [ot/(q — 1)], is given as in J5^| ) tuit/i 
go = g G [l,min{2,l + a}). Lei Xjv &e a symmetric 
Qa.k -independent (for some — oo < k < [a/(q— 1)] and 
a G (0,2])J random variables all having the same prob- 
ability density function fix) satisfying the conditions of 
Proposition \II.3\ 
Then the sequence 



Z N — 



X 



N 



is q a ,k- convergent to a iq a ,k—i, a) -stable distribution, as 
N -too. 

Proof. The case a = 2 coincides with Theorem 
1 of |7|. For k = 0, the first part of Theorem (q- 
convergence) is proved in Section 5 of the present pa- 
per. The same method can be applied for k ^ 1. For 
the readers convenience we proceed the proof of the first 
part also in the general case, namely for arbitrary k. 
Suppose that < a < 2. We evaluate F qa k (ZN). De- 
note Yj = Xj/s N iq a ,k),j = 1,2,..., where s N iq a<k ) = 

ifi qa k , a N) al2 - q '*.K\ Then Z N = Y 1 + ... + Y N . Again us- 
ing the relationship F q [aX]i£) = i^[X](a 2 "«£), we ob- 

■r). Further, it follows 

... and the associativ- 



tainF qak iY 1 ) = F qak [f]i- L — 

from q a ^.-independence of Xi,X2, 
ity property of the g-product that 



F ga , k [Z N m = ^ a F q .[f]i 



-), (61) 



' (Nck,aN) 

the right hand side of which exhibits the q Qj fe-product of 
N identical factors F Qa k [/]( 



,N)' 



Hence, making 



VI. SCALING RATE ANALYSIS 

In paper 7] the formula 

3 - qk-l 



Ik 



( 



Uq^- 1 -" 



2 -ik-i 



(64) 



was obtained for the g-Gaussian parameter f3 of the at- 
tractor. It follows from this formula that the scaling rate 
in the case a = 2 is 



5 = 



1 



qk-i 



Qk+l, 



(65) 



where qu~i is the g-index of the attractor. Moreover, if 
we insert the 'evolution parameter' t, then the translation 
of a g-Gaussian to a density in Q q [a] changes t to 1?l a . 
Hence, applying these two facts to the general case, < 
a < 2, and taking into account that the attractor index 
in our case is q* a k _ 1 , we obtain the formula for the scaling 
rate 



(66) 



«(2-C fe -i) 



In accordance with Proposition IV.21 2 — q* fe l = 
VCfe+i- Consequently, 

2 2 a-(fc-l)(g-l) 
d ~ a^k+i- aa _ (k + 1){q _ iy (67) 

Finally, in terms of Q = 2q — 1 the formula (|67l) takes 
the form 



2 2a-(fc-l)(Q-l) 



(68) 



a 2a- (k + l)(Q- 1) 

In Q it was noticed that the scaling rate in the non- 
linear Fokker-Planck equation can be derived from the 
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model corresponding to the case k = 1. Taking this fact 
into account we can conjecture that the scaling rate in the 
fractional generalization of the nonlinear Fokker-Planck 
equation is 



l-Q' 



which can be derived from (|68l) setting k = 1. In the case 
a = 2 we get the known result S — 2/(3 — Q) obtained in 



VII. ON ADDITIVE AND MULTIPLICATIVE 
DUALITIES 

In the nonextensive statistical mechanical literature, 
there are two transformations that appear quite fre- 
quently in various contexts. They are sometimes re- 
ferred to as dualities. The multiplicative duality is defined 
through 

fi(q) = l/q, (69) 

and the additive duality is defined through 

v{q) = 2-q. (70) 

They satisfy fi 2 = v 2 = 1, where 1 represents the iden- 
tity, i.e., l(q) = g,Vg. We also verify that 

(H m M m = M m (H m = 1 (m = 0, 1, 2, ...) . 

Consistently, we define (/i^)~ m = (^/i) m , and (vfi)~ m = 
Also, for m — 0, ±1, ±2, and Vg, 

= m-(m-l) g = q + mjl-q) 



v{piv) m (q) = - 
and 

(HX?) 



m + 1 — to q 1 + m(l — g) 



. + 2 - (to + 1) 5 2-g + m(l-g) 



to + 1 — mq 



1 + m(l — q) 



-m, + 1 + m q 1 — m(l — q) 



-to + (to + 1)? 9 — to(1 — g) 



We can easily verify, from equations ([Ml) and (|7T|) . 
that the sequences 52, n = 0, ±2, ±4, ...) and gi jn (n = 
0, ±1, ±2, ...) coincide with the sequence (^Lv) m (q) (m = 
0,±1,±,2,...). 



VIII. CLASSIFICATION OF (q, q)-STABLE 
DISTRIBUTIONS AND SOME CONJECTURES 



The g-CLT formulated in [7J states that the appro- 
priately scaling limit of sums of g^-independent random 
variables with a finite (2% — l)-variance is a ^-Gaussian, 



which is a g^-Fourier preimage of a (^-Gaussian. Here qk 
and ql are sequences defined as 

2g-fc(g-l) 



*(?-!) 



and 



qt = Qk-i, k = 0,±1, .... 
Schematically g-CLT in [7J can be represented as 

{/ : <T2 qk -i{f) < 00} A G qk [2] 3 &.[2]. 



(72) 



We have also noticed that g-CLT can be described by the 
triplet (P at t,Pcor,Pscl), where P att , P cor and P sc i repre- 
sent parameters of the attractor, the correlation and the 
scaling rate, respectively. We found that (see details in 
0) for g-CLT this triplet 



(Patt, Pcor, Psci) = (?fe-l j (Zfc+l)' 



(73) 



Schematically Theorem 3 of the current paper can be 
represented as 

C q [a] A £ g [a] A 0^ [2] , < a < 2, (74) 

where is the set of (g, a)-stable distributions, Q„l [2] 
is the set of g L -Gaussians with index q L defined as 



Q L =qa(q) = 



3 + (2g- l)a 



1 + a 



Recall that the case a = 2 was peculiar and we agree to 
refer to the scheme ([72]) in this case. 

Theorem 4 generalizes the g-CLT (which corresponds 
to a = 2) to the whole range < a < 2. Schematically 
this theorem can be represented as 



1a,k 



g qlk [2],0<a<2, (75) 



generalizing the scheme (|72p . The sequences q a .k and 
g* fe in this case read 

gg + fc(l-g) 
9a, fc = , "v, k = 0, ±1, 



and 



9a, fc 



1 - 



c* + fc(l -g) 
2(1 -g) 



a + ft(l -g) 



fc = 0,±1, 



Note that the triplet (P a tt, Pcor, Psci) mentioned above, 
in this case, takes the form 

{Patt, Pcor, Psci) = {q* a ,k-l, Qa,k, (2/ a )5a,fc+l)j (76) 

recovering the triplet (|73p in the case a = 2. 

In connection with the above discussion about triplets, 
we note that the existence of a g-triplet, namely 
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(<lsen, Qrei, Qstat) , related respectively to sensitivity to the 
initial conditions, relaxation, and stationary state was 
conjectured in [36]. Later it was observed in the solar 
wind at the distant heliosphere [13, [38[ . The triplet in 
([73)) obtained theoretically might be useful hint for its 
understanding. 

Finally, unifying the schemes ([71)1 and ([75)1 we obtain 
the general picture for the description of (q, a)-stable dis- 
tributions: 



(77) 



where 



la.k = 3a(<2a,fe) 



3 + (2g a , fc - l)a 
1 + a 



finite i::lassjc vananc e 



1.4 - 

1.2- 



FHK^ I Q-Cauchy li 




7 WW 



FIG. 1: (Q, a)-rcgions. 



In Fig. Q]the dependence of q L and q* on parameters 
(Q, a) £ Q in the case k = is represented. If Q — 1 and 
a = 2 (the blue box in the figure) , then the random vari- 
ables are independent in the usual sense and have finite 
variance. The standard CLT applies, and the attractors 
are classic Gaussians. 

If Q belongs to the interval (1,3) and a = 2 (the blue 
straight line on the top), the random variables are not 



independent. If the random variables have a finite Q- 
variance, then g-CLT [7| applies, and the attractors be- 
long to the family of g*-Gaussians. Note that q* runs 
in [1,5/3). Thus, in this case, attractors (<?*-Gaussians) 
have finite classic variance (i.e., 1- variance) in addition 
to finite g*-variance. 

If Q = 1 and < a < 2 (the vertical green line in 
the figure), we have the classic Levy distributions, and 
random variables are independent, and have infinite vari- 
ance. Their scaling limits-attractors belong to the family 
of a-stable Levy distributions. It follows from (j44[) that 
in terms of g-Gaussians classic symmetric a-stable distri- 
butions correspond to ^§/3<qo,Gq @. 

If < a < 2, and Q belong to the interval (1,3) we 
observe the rich variety of possibilities of (q, a)-stable 
distributions. In this case random variables are not inde- 
pendent, have infinite variance and infinite Q-variance. 
The rectangle {1 < Q < 3; < a < 2}, at the right 
of the classic Levy line, is covered by non-intersecting 
curves 



C qI .={(Q,a) 



3 + Qa 
a + 1 



q L }, 5/3 < q L < 3. 



This family of curves describes all (Q, a)-stable distribu- 
tions based on the mapping (|74l) with g-Fourier trans- 
form. The constant q L is the index of the g L -Gaussian 
attractor corresponding to the points (Q, a) on the curve 
C q L . For example, the green curve corresponding to 
q L = 2 describes all Q-Cauchy distributions, recover- 
ing the classic Cauchy-Poisson distribution if a — 1 (the 
green box in the figure). Every point (Q, a) lying on the 
brown curve corresponds to q L — 2.5. 

The second classification of (Q, a)-stable distributions 
presented in the current paper, and based on the mapping 
([75[) with q* -Fourier transform leads to a covering of Q 
by curves distinct from C q L . Namely, in this case we have 
the following family of straight lines 



{(Q,a) 



4a 



Q + 2a 



— = 3- 9 *}, l<q* <3, (78) 



which are obtained from (|56l) replacing n = — 1 and 
2q — 1 = Q. For instance, every (Q, a) on the line F-I 
(the blue diagonal of the rectangle in the figure) identi- 
fies q*-Gaussians with q* = 5/3. This line is the frontier 
of points (Q,a) with finite and infinite classic variances. 
Namely, all (Q,a) above the line F-I identify attractors 
with finite variance, and points on this line and below 
identify attractors with infinite classic variance. Two 
bottom lines in Fig. [TJreflect the sets of q* corresponding 
to lines {1 < Q < 3; a = 2)} (the top boundary of the 
rectangle in the figure) and {1 < Q < 3; a = 0.6)} (the 
brown horizontal line in the figure). 

Some conjectures. Both classifications of (Q, te- 
stable distributions arc restricted to the region Q = { 1 < 
Q < 3,0 < a < 2}. This limitation is caused by the 
tool used for these representations, namely, Q-Fourier 
transform is defined for Q > 1. However, at least two 
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facts, the positivity of // g>a in Proposition 111.31 for q > 
max{0, 1 — 1/a} (or, the same, Q > max{— 1, 1 — 2/a}) 
and continuous extensions of curves in the family C q r , 
strongly indicate to following conjectures, regarding the 
region {Q < 1} on the left to the vertical green line (the 
classic Levy line) in Fig. [TJ In this region we see three 
frontier lines, F-II, F-III and F-IV. 

Conjecture VIII. 1 The line F-II splits the regions 
where the random variables have finite and infinite Q- 
variances. More precisely, the random variables corre- 
sponding to (Q, a) on and above the line F-II have a fi- 
nite Q-variance, and, consequently, q-CLT Jj/ applies. 
Moreover, as seen in the figure, the q L -attr actors corre- 
sponding to the points on the line F-II are the classic 
Gaussians, because q L = 1 for these (Q,a). It follows 
from this fact, that q L -Gaussians corresponding to points 
above F-II have compact support (the blue region in the 
figure), and q L -Gaussians corresponding to points on this 
line and below have infinite support. 

Conjecture VIII. 2 The line F-III splits the points 
(Q,ct) whose q L -attr actors have finite or infinite classic 
variances. More precisely, the points (Q, a) above this 
line identify attractors (in terms of q L - Gaussians) with 
finite classic variance, and the points on this line and 
below identify attractors with infinite classic variance. 



Conjecture VIII. 3 The frontier line F-IV with the 
equation Q + 2a — 1 = and joining the points (1,0) and 
(—1,1) is related to attractors in terms of q* -Gaussians. 
It follows from |7£| ) that for (Q,a) lying on the line F- 
IV, the index q* = — oo. Thus the horizontal lines cor- 
responding to a < 1 can be continued only up to the line 
F-IV with q* 6 (-co, 3 — q + 2° _i ) ( see the dashed hori- 
zontal brown line in the figure). If a — > 0, the Q-interval 
becomes narrower, but q* -interval becomes larger tending 
to (— co, 3). 

Results confirming or refuting any of these conjec- 
tures would be an essential contribution to deeper un- 
derstanding of the nature of (Q, a)-stable distributions, 
and nonextensive statistical mechanics, in particular. 

Finally we note that Fig. [1] corresponds to the case 
k = in the description (|77|) . The cases k ^ can be 
treated in the same way. 
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